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CHORD WISE AND COMPRESSIBILITY CORRECTIONS TO SLENDER- WING THEORY 1 


By Harvard Lomax and Loma Sluder 


SUMMARY 

Corrections to the solutions given by slender-winy theory jar 
the lift distribution on triangular and rectangular wings of low 
aspect ratio are obtained in two steps. First, slender wing 
theory is used to find the load distribution over a wing of given 
shape. Second, the spanwise variation of the loading so ob- 
tained is left unchanged but the chordwise variation is modified 
by satisfying an appropriate integral equation. Results are 
shown for flat-plate wings and, in the case of the subsonic, 
triangular wing, a comparison is made with other theoretical 
solutions and experimental results. 

INTRODUCTION 

The calculation of loading on three-dimensional lifting 
surfaces is a fundamental problem in aerodynamic research. 
The complexity of the problem has led to the development 
of certain simplified theories by means of which the loading 
on special types of plan forms can be estimated quickly. 
The amount of error which these estimates contain is of 
considerable interest, as are methods which will tend to 
correct such errors without undue labor. 

Slender-wing theory applies to one, such simplified body 
of analysis. There are two basic assumptions of this theory. 
One, the angle of attack is small enough so that the vortex 
sheet lies in the plane of the wing and the boundary condi- 
tions for the wing can be projected onto a horizontal plane 
parallel to the direction of the free stream; and the other, 
that either the chordwise gradient of velocity is small enough 
or the free-stream Mach number is close enough to unity 
that the linearized partial differential equation which governs 
the fluid flow becomes Laplace’s equation in a plane trans- 
verse to the free-stream direction. References 1 through 6 
are examples of papers developing slender-wing theory. 

The magnitude of the error of such a theory, in the case 
of subsonic flow, is indicated by observing solutions for 
triangular wings. Slender-wing theory gives a finite value 
for the loading along the trailing edge. Proper inclusion of 
the chordwise and compressibility effects results in solutions 
that satisfy the Kutta condition and make the loading fall 
to zero at the trailing edge. It is the purpose of this report 
to study such modifications. 

The corrections due to the chordwise and compressibility 
effects are obtained in the following manner: First, an in- 
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tegral equation is set up relating the shape of the wing 
surface to the lift distribution; second, this integral equation 
is solved for the given shape under the assumption that the 
chordwise velocity gradients are small or that the free-stream 
Mach number is unity; and finally, the integral equation is 
reinspected for the same wing shape, this time with the 
spanwise lift distribution fixed at the variation just obtained 
and with the chordwise variation as the unknown and the 
Mach number terms included. 

Results are presented and discussed both for triangular 
and rectangular, flat-plate plan forms in both subsonic and 
supersonic flow. 
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LIST OF IMPORTANT SYMBOLS 


aspect ratio (jr) 


span of wing measured normal to plane of 
symmetry 
root chord of wing 

lift coefficient 

pitching-moment coefficient 

''pitching moment about leading edge or apex' 


( E 


aSc a 


) 


incomplete elliptic integral of the second kind 
with, argument t and modulus k 


c^rvw*] 


complete elliptic integral of second kind, E (1 ,£) 
incomplete .elliptic integral of the first kind 
with argument i and modulus k 


f {t,k) JoVd-oa-w)] 


complete elliptic integral of first kind, F{l,k) 
for triangular wing, slope of leading edge relative 
to plane of symmetry 
free-stream Mach number 

( Y* } 

\speed of sound in free stream/ 
static pressure 
Pi—P* 

free-stream dynamic pressure ^ PoFo*) 

semispan of rectangular wing 
area' of wing 
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do 

Po 

a, <r„ 
<P 

Aip 
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u 


perturbation velocity component in the direction 
of the x axis 

U u —Ui 

free-stream velocity 

perturbation velocity component in the direction 
of the 2 axis 
-Vo a 

Cartesian coordinates of an arbitrary point 
Cartesian coordinates of source or doublet 
position 

Co 

angle of attack 

vri=22?r 

/3m 

density in free stream 
doublet weighting factors 
perturbation velocity potential 
<Pu~ <Pl 

SUBSCRIPTS 

conditions on lower surface of wing (at 2=0—) 
conditions on upper surface of wing (at 2=0+) 


In linearized theory this jump can be related to the loading 
coefficient A p/g by the equation 




( 2 ) 


Further, the velocity potential <p can be found by the relation 



Operating on equation (1) in this manner and interchanging 
the order of integration gives 
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(3) 


which represents, physically, a distribution of elementary 
horseshoe vortices. 

The effect of compressibility in a linearized study of lifting- 
surface theory can only enter through the use of /3. Setting 


THE INTEGRAL EQUATIONS 

SUBSONIC 


<r=l + 

v 


? (* 


X — Xj 

— *0*+ PHy — yi)*+ p V 


(4) 


Triangular plan form. — A general solution of Laplace’s 
equation which is suited to problems in linearized subsonic 
wing theory (given, e. g., in reference 7) is that which relates 
a velocity potential or perturbation velocity to the value of 
its jump across a given surface. For the lifting triangular 
wing shown in figure 1 this can be written 

p -k&i— m 

u 4 x Jo w 


it is seen that a is the only term in equation (3) which con- 
tains j8. This term has an interesting interpretation in the 
light of the study which has been made at sonic speeds. At 
M a =l (i. e., /9=0), o- has either the value 2 or 0, depending 
on whether aq is less or greater than x. Hence, for il/ 0 = 1, 
equation (3) becomes 


4x Jo J-mxi(y— 




ytf+z* 


(5) 


where /3=V|1 — M*\, u is tho perturbation velocity parallel 
to the x axis, and Au is the jump in u over the wing plan form. 


Now reversing the order of integration and using the define- 
tion implied by equation (2), namely, 



t'inuRS 1.— Trlaagular-wlng ooordinate system. 


Ap_ 2 dA«> 
g ~ Vo dr, 

gives finally 

_ g f "* Apdyi 

* Zrj-mxiy—yiY+z* 


( 6 ) 


Equation (6) has been studied in reference 5 as the funda- 
mental equation for slender wings or wings flying at near 
sonic speeds. It is an equation which gives the solution for 
the velocity potential in a tlLree-dimensional flow in terms of 
two-dimensional doublets, the two dimensions being at right- 
angles to the free-stream direction. A solution of such a 
nature is immediately implied by the physical character of 
both sonic wing theory, in which the .Mach cone has degener- 
ated to a Mach plane, and slcnder-wing theory, in which the 
wing is so slender that the chordw'ise gradient of velocities 
can bo neglected compared to the vertical and lateral 
gradients. 

By comparison of equation (5) with equation (3), it is 
seen that the term a can be interpreted as a factor which 
corrects the slendcr-wing-theory results as given by equation 
(5) for the effects of chordwise gradients in velocity and 
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compressibility. By consideration of the effect at one point 
of the distribution of doublets over the wing, this correction 
can be visualized as a reweighting of the two-dimensional 
doublets according to their position relative to the point. 
Figure 2 indicates the variation of a across the span at various 
chord stations for (3=0.6. Observe that the doublets ahead 
of the point at which the potential is to be determined are 
still weighted far more heavily than those behind the point. 
The effect of considering (3 different from zero, however, is 
to reduce the extreme difference in weight occasioned at 
£=0 so that the doublets be hin d a given point do have some 
effect on the induced velocities there, and the doublets ahead 
of a point induce a somewhat smaller disturbance than 
before. Since the strength of these weighted two-dimensional 
doublets is given by the magnitude of the three-dimensional 
loading, their strength is zero everywhere off the wing plan 
form including the area behind the wing occupied by the 
vortex wake. 

Two different methods for the further reduction of equation 
(3) will be considered. The first method involves finding the 
vertical induced velocity for points along the x axis, while 
the second involves finding the average vertical induced 
velocity along the span at a given chord station. The first 
method must be discarded for triangular wings because of 
difficulties around the apex; the second, however, proves to 
be satisfactory. The simplification obtained by considering 
the vertical induced velocity for points along the x axis will 
be considered later in connection with the rectangular wing. 

Since it is easier to consider first the averaging process, 
the operator 


lim _j_ A 
*-<! 2 rax dz 



- r “fv r i- ( v ° rdx r dyi 

w 8?Jo dXl J 8xJo dXl }-^2mx (x-x.) 

r y’lar-rQ^iS^mr-yi) 1 VQ— x l ) t +(S i (mx+ yQ 1 ! „ 
[_ mx—yi ‘ mx-\-yi J ' ' 


where w is the average value of the vertical induced velocity 
along a given span. _ . . 

The solution for A p/y obtained from slender-wing theory 
can be written * 


A p 
2 


4rcom a r‘i 

F 0 -v / mhr7— y? 



( 8 ) 


where in that theory /i(xi/c 0 )=l. If the value of A pfq 
given by equation (8) is placed in equation (7), the resulting 
integral equation can be written in a simplified form if it is 
noted that 



'(x— r t )*+/3 2 (/nx — ytf 
mx—yi 


Jjx—XiY+PXmx+yiY 

mx+yi 


\mdyi 

■JirPxf—y' 


“ m J. 


mCt+r,) 


■^(x—XiY+p^dri 
tl\ — m^x *— — ij 1 


(9) 


where for the first term in the brackets the transformation 
H—mx—yi was used and for the second the transformation 
ri—mx+yi. Hence, equation (7) finally reduces to the 
following 



The solution* of equation (10) will be discussed in a later 
section devoted to triangular wings. 

Rectangular plan form. — If the plan form of the wing is 
rectangular as shown in figure 3, then equation (3) is modified 
slightly to the form 




{y-yif+z 1 


( 11 ) 


It is possible in this case to study the vertical induced veloc- 
ity for points along the x axis; that is, to find 2 )tpjbz by 
equation (11) and then set both y and z equal to zero. In 
order to do this a special notation is employed. Thus, if 
the indefinite integral of/(y)/y* can be written (where /(y) is 
bounded at y=0) 

dy=J(y)+C 

then, by definition, 

£f>dy^(«W(-«) (12) 

By means of this definition it can be shown (see reference 8) 
that 


is applied to the weighted doublets, trz/[(y—y i) i +3 2 ], of equa- | 
tion 13) with the result that t 


> This solution follows from an analysis of equation (J) . See reference 5. 
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Fkidri 3.— Rectangn tar-wing coordinate system.-. 
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If 


A p ^ w 0 

q VqCq 



(14) 


integration by parts gives 


and equation (13) can be written 


~Vi 

~+Wy? 


(15) 



_h_ 

Z — Xi 



dx i 


(16) 


This integral equation has been derived previously by K. 
Wieghardt (reference 9) with regard to the rectangular-ving 
problem. The solution of equation (16) will be discussed 
in a later section devoted to rectangular wings. 


SUPERSONIC 

Triangular plan form. — In passing from subsonic to 
supersonic theory, we pass from the elliptic to the hyperbolic 
partial differential equation and in particular from Laplace’s 


« Th« solution for the rectangular wing given fcy sleader-wing theory Is that the load be 
zero across every spanwise strip a/t of the leading edge. To find the chordwise correction to 
suoh a theory, therefore, a spanwise distribution must be assumed. Since, however, slender- 
wing theory also requires an elliptical span loading for the boundary conditions of a rectan- 
gular wing to be satisfied, a reasonable choice Is that given by equation (14). 


equation to the w r ave equation. The solution which relates 
the perturbation velocity u at any point in the field to the 
loading on the wing can again be written in terms of an 
elementary horseshoe vortex distribution over the wing plan 
form. As in reference 8, this becomes 

. ff (z — Zi)Audx\dyi 

2t dx J J, [(y —yif-r z*\ y/(x — Xi) s — 0*(y — yi) s — /9*3* 

. . , Ap 2A« 

and since u=s— and -=-= 7 -> 
ox q Vo 

Voz rr ~ f dx ldyi r x-h 

^ 4 t J J t (y—y 0 a +2 s [_V(z— 2i)*— P*(y— vO*— 

where r is the area on the wing bounded by the edges and 
the trace of the Mach foreconc from the point x, y, z. Again 
the effect of compressibility appears only in the term within 
the brackets. Hence, setting 



V(x— Zi)*— 0*(y— y,)*— /3 s z ! 

o-p contains all of the Mach number effects at supersonic 
speeds. At M 0 =l, o-, = l, and since, by the definition of 
r, x t <x, it follows that at sonic speeds equation (17) also 
reduces to equation (5). Therefore the doublet distribu- 
tions represented by equations (17) and (3) are consistent 
at the speed of sound. 

In order that an exact parallel can be provided with the 
subsonic solution to the triangular wing, the average vertical 
induced velocity for points along a given span is again 
considered. It can be shown (reference 8) that in the plane 
of the wing 


v)—- 
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Vo 


■ 4 


4tt 


(*x r^i 

R.P. dx i dyi 


mx i 


(*— *>) Y 


(y — y i) V(z — xi)*— /S’(y — y ,)* 


(19) 


where the order of integration must be carried out as indi- 
cated (i. e., the integration with respect to yi must be made 
first). The letters I?. P. mean that the real part of the term 
is to be taken. Such a device can be used since the double 
integral must always be a pure real quantity in the area r 
(Ap/y is real everywhere on the plan form) and a pure 
imaginary quantity over the rest of the area indicated by 
the limits on the integrals (sec fig. 4). The average vertical 
induced velocity along the span may be obtained by applying 
to equation (19) the operator 


and since 


Ap 


2^fZ dV 

4«> 0 m s xj 




q Voi/mW-vS 

equation (19) may be written in the form 

— jKPmwojtjx) . Wq f J Zi/i/t(gi)dzt 
2 t "rxjo x—xi 


( 20 ) 


( 21 ) 


X—Xi 
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Fnura* 4.— Regions for which tho double Integral term In equation (I© b a pure real or pure 

Imaginary. 



The derivation of /, is similar to that used for equation (9) 
and leads to the expression 

I a =m R.P. lY-fv'dri ' (22) 

.Jmd-i,) 7jV[?n(zi+x)— t;l ar)+77] 

It is possible to find an exact solution for /*(x) by means 
of equation (21), but the discussion of this analysis is reserved 
for a subsequent section. 

Eectangular plan form. — Equation (19) can also be used in 
the case of a rectangular plan form by an appropriate change 
in limits; thus, 



where again it should be stressed that the order of 
integration cannot be reversed. The regions in which 
the integration yields real or imaginary results are 
shown in figure 5. As in the case of the subsonic rectangular 
wing, the value of to will be obtained only along the x axis 
so that y in equation (23) can be set equal to zero. The 
loading will be assumed to have a form 





(24) 


which is similar to that used in the subsonic case except that 
the reference length is now the semispan instead of the chord. 
Such a difference is reasonable since in the supersonic case 
the position of the trailing edge cannot affect the loading 
on the wing. 

Finally, therefore, when y=0 equation (23) becomes 
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Pure v 

imaginary 


Fioubi 5. — Regions for which the double Integral term fn equation (») fs pure real or pure 

imaginary. 


(25) 


(26) 

The solution of equation (25) is deferred to a subsequent 
section. 

LOADING ON WINGS 

The previous section was devoted to the development of 
the integral equations which are to be studied for the two 
types of plan forms in subsonic and supersonic flight. In 
order that this study can proceed in a natural manner, the 
arrangement of the presentation has been changed so that 
the plan form is the principal division and the speed is 
subsidiary. 

TBIA.NGULA.R WIXGS 

Supersonic case. — The decision to solve for the loading 
on the supersonic, triangular, flat plate by analyzing equa- 
tion (21) was not an obvious one since the exact solution of 
the linearized partial differential equation for this case has 
already been obtained. (See, e. g., references 10, 11, and 
12.) Thus it is known before starting that the value of 
/g(x) in equation (21) must be 1/E where E is the complete 
elliptic integral of the second kind with modulus ^1—771*0*. 
However, these solutions were obtained by an entirely 
different procedure so that by solving equation (21) and 
comparing the two results a check on the accuracy of the 
method is obtained. Furthermore, when the subsonic 


\80j Jr8j3j 0 X — Xi 


where is given by the equation 


h=R.P.J\ 


/(x-x i)*— 18*1/1* 
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problem is analyzed the same general procedure will be 
followed and the results can then be accepted with greater 
confidence. 

The fust step in the solution of equation (21), in which w 
has been set equal to w a since the wing is a flat plate, is to 
change variables by the transformation £ l =x l /'x. Tliis gives 

i f i m (27) 


in the equation for I t , the transformation m.r} X = -qjx was used 
so that 

r . f 1+ft V(l-fr)«-g»CTW . ' 

J 1_ fi T)1 V(h+ 1 — » 7 i) (fl — 1 + 171) 

which is completely independent of x. The partial deriva- 
tive of both sides of equation (27) with respect to x gives 


/.'(*) 


2 p feV/Cgfelrffc 

ir 2 j3m Jo l—|i 


h&) 


(28) 


‘Equation (ill) is a special form of Abel's integral equation, 
the unique inversion 4 of which is, in this case, /,(x) = l, 
This is easily verified by direct substitution. 

The simplicity of this result is not accidental, of course, 
since the valuo of / s ( x) was originally introduced by equa- 
tion (20) as a correction factor to the slender-wing-theory 
solution. 

Subsonic case.— The study of the triangular wing pre- 
sented in the preceding section was made first at arbitrary 
supersonic Mach numbers and then at a Mach number equal 
to 1. In keeping with this order of decreasing speed, the 
subsonic flat plate will be studied first at sortie speed and 
then for general subsonic Mach numbers. 

An inspection of equations (9) and (22) is sufficient to 
show that (Ii)fl*_o is equal to (/*)ffm-o- Hence, equation 
(30) can be substituted into equation (10) and there results 
(since again w is set equal to w 0 ) 


1 

p © ix> , 


2 x 2 ; 

LJ v * 2 -*! 2 "*■ 

Jo aV 


Equation (28) is a homogeneous linear integral equation. 
The solution to equation (28) is simply ji(x)—Q or, what is 
equivalent, /j(z) equals a constant, (/a)o say. By means of 
equation (27), this constant can be evaluated. Hence, 



which represents the solution to the problem. The integral 
7, was calculated analytically as in appendix A, and then 
the value of (/ 3 )o, as given by equation (29), was determined 
by numerical integration. For )3m=0.8 the result of this 
computation was 0.708; wheieas the true value given by 
l/E is 0.705. 

Equation (27) can also be solved when the wing is slender 
with respect to the Mach cone by considering fim to be 
small. Setting fim=Q yields 

( s)»m 0 ( l)iV(£i+l — 111) (Si— 1+171) 

artd this is readily evaluated to give 

( h)fn - 0= T===f ( 30 ) 


The integral equation reduces to 

1== r 1 6/a(gfe)<*fe 

Jo 1 — £i ! 

which by a re transformation of variables Xi=x£i becomes 
x= C* xJi(3:i)dx i 

Jo ^X t ~Xi i 


and this reduces immediately to 



It is now obvious that equation (31), which was derived 
from supersonic wing theory, and equation (32), which was 
derived from subsonic wing theory, are identical. Clearly 
this establishes the continuity of the theory in passing from 
the supersonic to the subsonic regimes. 

The study of the general subsonic case leads eventually 
to the numerical solution of an integral equation. However, 
an idea of the qualitative form which this solution must 
assume can be gained by some preliminary analysis. 

First write equation (10) in the form 



f 10 Tr z ij_i( x t)dxi | f l X ^ 1 


_J V^O*— Xl* J U 

10 — 3-2 


(33) 


where Xq=x/c 0 , x i =Xijco, and w/w a = 1. The evaluation of 7i 
is given in appendix A, and a plot of £/i against Xi/x 0 for 
(/3m) J equal to 0, 0.05, 0.10, and 0.20 is shown in figure 0. 
Obviously equation (33) is a singular integral equation. 
We have already seen that for /Sm=0 it is an Abel type 
integral equation with a K-power singularity at the upper 


* If Abel's equation is written In the form 


Its Inversion Is 




■ and this inversion Is unique for functions g(x) that fulfill the condition 
Urn 

e-+Q -0 . a£x <f (a) 


( 31 ) 


Since the solutions to r the velocity throughout the flow field must also satisfy condition CD. 
the value of g(x) given above la unique In the class of functions available. 
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limit. For finite values of 0m, however, I t is seen from 
figure 6 to be bounded and nonzero in the interval 
00 ; hence, for such 0m, equation (33) is a combi- 
nation of two types of singular integral equations, the Abel 
type and the Cauchy type, the latter having a first-power 
singularity in the interval of integration. 

Experience with the Cauchy type integral equation 


Jo y— 


g(x)dx 
x 


( 34 ) 


winch arises, for example, in the study of subsonic lifting- 
line theory and two-dimensional airfoil theory is useful in 
the present problem. Thus, the solution to equation (34) 
is not unique (even when g(x) is restricted according to 
equation (a) in footnote 4) unless some additional condition 
is given. Such a condition might be the requirement that 
tf(l)=0; this would correspond in two-dimensional airfoil 
theory to the specification of the Kutta condition at the 
trailing edge. Further, it is known that if the condition 
,?(1)«=0 is satisfied in equation (34), then g(x) tends to 
infinity as x approaches zero. 

Since the solution to equation (33) must also satisfy the 
Ivutta condition, the above discussion leads one to anticipate 
for the shape offi(x<>) a curve something like that shown in 
figure 7 (a). Ou the basis of such qualitative knowledge, a 
simple numerical procedure was set up and used to calculate 
the solution to equation (33). This procedure, which is 


based on the assumption that/i(r 0 ) is constant over each of 
nine equally spaced intervals, is presented in appendix B. 
The results of the analysis for (0m) 1 equal to 0.1 (i. e., 
jSA = 1.26) are shown in figure 7 (b). 

In order to check the results derived from the method just 
mentioned, equation (33) has been solved in an alternative 
manner. In this second ^approach it is assumed that /i(^o) 
can be approximated by a second degree polynomial which is 
multiplied by two factors, one that vanishes for r 0 =l and the 
other that tends to infini ty as x 0 vanishes. 

First consider the behavior of / t (x 0 ) for small x Q . A brief 
s.tudy of the solutions to the two integral equations repre- 
sented by equations (32) and (34) will illustrate how f L (xo) 
can be analyzed in the vicinity” of the origin. First consider 
the Gauchy type integral equation (34) and assume that for x 
small g(x) can be expressed in the form 


•c »— o 


1>5>0 

Then equation (34) becomes 

~ f 1 dx - C'-x'-'dx 

“Jo x\y—x) jS a *J Q x—y 

which in the limit as y —> 0 reduces to 



Fwrsi 6. — Variation of /t/2 with 
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Pinmts 7.— Variation of /i'w'lth zt. 


or, setting x—y%, to 


UmLC m 

»-° 2/*Jo 
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The term on the left is indeterminate and implies 


lim C 

*-°Jo 


Uv 


€*(1-0 



dX 


cot [(1— a)ir] = 0 


This equality is satisfied, since 5 is to be greater than zero and 
less than ], only by the value 5=}i. The fact that this must 
be the exponent of 1/x in the expression for g(x) can be verified 
by inspecting the known solution to equation (34). 

A similar analysis applied to Abel’s integral equation in the 
form given by equation (32) shows for that case J must be 
zero which, again, agrees with the known inversion. 

Finally, in appendix C this same approach is used to dis- 
cover the initial behavior ofji(zo) in equation (33). Figure 8 
presents the results of this analysis throughout the range of 
0A for which I t was calculated. 

As was already mentioned, the variation of /i(x 0 ) ha tlio 
vicinity of x 0 =l, that is, near the trailing edge of the wing, 
is fixed by the Kutta condition. A useful statement of. this 
condition that holds in both subsonic and supersonic flow is 



PA 


TTfltms 8.— Variation of t with reduced aspect ratio PA. 

that the magnitude of the loading at the trailing edge must 
not be. infinite. However, the only pertinent solution to 
equation (34), and hence to equation (33), that is not infinite 
at x=l is that which is identically zero there. Further, ns 
the trailing edge of the triangular wing is more and more 
closely approached, it is reasonable to .expect that the shape 
of the load distribution in its vicinity approaches that in the 
vicinity of a simple two-dimensional wing trailing edge, Lite ef- 
fect of the wing plan form vanishing as the ratio of the dis- 
tance from the leading to trailing edge tends to infinity. On 
the basis of these considerations, /i(x 0 ) should approach zero 
as the term (1 — Xo) w approaches zero when x a tends to 1 . 

It is also apparent, however, that/i(xo) equals 1 for all x a 
between zero and 1 when the Mach number is unity. Fur- 
ther, as M 0 approaches 1 or as the angle of sweep approaches 
90°, the effect of the plan form on the shape of the loading 
near the trailing edge becomes increasingly important or, in 
other words, the trailing edge of the triangular wing must be 
more and more closely approached before the shape of the 
two-dimensional load distribution is simulated. An ex- 
ponent to (1 — Xo) which satisfies this requirement as well as 
those in the preceding paragraph is 7/2(7 -(-1— £0) where 7 is a 
function of Maeh number and vanishes as M 0 — >4. 

Finally, therefore, it is assumed that/i(x 0 ) can be expressed 
by the equation 

fl— 

/i(zo)=- 2 ~~r (a 0 +a,x 0 +ajx 0 *) (35) 

Xo 

The values of the constants a 9 , a 1( c tj, and 7 are determined by 
satisfying equation (33) at four chordwise stations and are 
given for (/3m) * equal to 0, 0.1, and 0.2 in table I. The 
accuracy to which those solutions satisfy the integral equation 
is indicated by table II. 

A comparison, for (0m) a =O.l, between the solution given 
by equation (35) and that derived by the method outlined in 
appendix B is shown in figure 7 (b). The strictly numerical 
method presented in the appendix was used for three different 
interval spacings, results for which are indicated in the fig- 
ure by the symbols. Presumably the accuracy of the method 
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increases with increasing number of intervals used. A study 
of the figure shows that the numerical method is apparently 
approaching 5 the solution given by equation (35). Subse- 
quent values in this report are based on the solutions repre- 
sented by the latter equation. 

The final curves for/ifo) are shown in figure 9. A discus- 
sion of the integrated values of the loading will be given later. 

RECTANGULAR WINGS 

The discussion of the triangular wing was divided accord- 
ing to the Mach number. The same division will be used in 
tliis section, starting with the discussion of the results for 
supersonic speeds, then with that for both supersonic and 
subsonic theories at sonic speeds, and closing with a discus- 
sion of the subsonic development. 



Supersonic case. — The solution of equation (25) will give 
the loading on a rectangular wing flying at a supersonic 
Mach number. The evaluation of the integral 7 t is carried 
out in appendix A where it is shown that /* can be expressed 
in terms of complete elliptic integrals of the first and second 
kinds. Having the expression for J 4 , a numerical solution 
may be obtained for (see appendix B). Figure 10 

shows a plot of (a factor representative of the chord lift 
distribution) for a flat plate wing as a function of xfsl 3, 
the ratio of the distance back from the leading edge to the 
magnitude of the reduced semispan. The value of / 4 given 
by equation (25) can be checked in the interval 0 < (x/s/S) <2 
because the exact solution to the complete linearized partial 
differential equation can be readily obtained there. The 
comparison is given in figure 10. The fairly rough agreement 
shown is not surprising since equation (25) is derived on the 
assumption that the spanwise distribution of load is elliptical 
at every chord station, and certainly tliis assumption is 
least accurate in the interval, where the comparison with 
the exact results is made. The area under the exact and 
approximate curves in figure 10, between the initial value 
and that at which / 4 =0, is nearly the same. (See the next 

* The maximum difference In total lilt is already leas than 8 percent (based on the lover 
value) and the difference tn center of pressure Is negligible . 


section on Aerodynamic Characteristics.) The integrated 
value of/ 4 as given by equation (25), therefore, can be used 
for x/s/9> 2. 

As for the qualitative nature of the variation, figure 10 
shows that the loading on a narrow rectangular wing flying 
at supersonic speeds falls linearly to zero, becomes negative, 
and then oscillates between negative and positive values, 
the amplitude of the oscillation being so heavily damped 
that after the third change in sign the magnitude is practically 
zero. 

It should be noticed in studying the results of figure 10 
that the entire resultant lift of the wing is concentrated 
in the interval 0<(x/«/3)<2. But as the Mach number 
approaches 1 this interval approaches zero, and the entire 
lift of the wing is carried in a strip along the leading edge.® 
Such a solution violates, in the vicinity of the leading edge, 
the assumption on which the theory is based and should 
be considered only as a theoretical limit. 

Results for the lift and pitching moment on the rectangular 
wing will he developed in a later section. 

Subsonic case. — The study of the subsonic rectangular 
wing stems from equation (16). The first step in the analysis 
of the equation will be to consider its solution at fis—O 
and show that this is continous with the supersonic results 
there. 

The value of h can be written (equation (15)) as 

(x-xO’+flVy* , 

1=7 d v 

and for ps= 0 this becomes 

l (X — Xi)v Xi<X 

(— (x — Xi)t Xi>x 

and hence equation (16) can be written 

— 5JX3)* 1 (36) 



* This result also folloivs by Inspecting equation (25) for the values £l»0. 



272483—54 68 
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Equation (36) is identical with the form of the supersonic 
equation (25) at )3s= 0 so that once again the continuity 
of the subsonic and supersonic theories at the sonic speed 
range is established. Furthermore, equation (36) shows that 
if w/w a is constant then Ji(x 0 ) must be zero everywhere 
except at points where it can be represented by a pulse, 
the integral of which has a finite magnitude. From the 
supersonic discussion, it is clear that one such pulse exists 
and is located at the leading edge. 

The evaluation of h for /3«>0 is given in appendix A. 
The numerical solution to equation (16), assuming the Kutta 
condition at the trailing edge, is given in appendix B for 
values of reduced aspect ratio pA equal to 0.33, 1.0, 1.5, 
and 2.0 For an aspect ratio equal to 2, these values corre- 
spond to Mach numbers of 0.986, 0.866, 0.662, and 0, respec- 
tively. The results of the computations are shown in figure 
1 1 where the chordwisc lift distribution factor/ s (zo) is plotted 
against x-o for the various values of PA. By comparison of 
figure 11 with figure 10, it can be seen that in the subsonic 
case the loading drops monotonicallv from infinity at the 
hading edge to zero at the trailing edge and does not oscil- 
late in the afterportion, as in the case of the supersonic wing. 

When p equals one, these results can be compared with 
those obtained by Wieghardt and presented in reference 9. 
Figure. 12 shows the comparison for two values of the aspect 
ratio. Curves are also shown in - the figure for the loading 
obtained by using the method given in appendix B but by 
satisfying the integral equation at only six and three points. 
The latter curve is in better agreement with Wieghardt’s 
result and, since Wieghardt (although using a different 
method involving Birnbaum functions) used only four points, 
tills may account for the discrepancy between the final 
results of this report and those of Wieghardt. 



Fiqobje II. — Variation of cbordwlse correction factor f\ for subsonic rectangular wing. 



Fiuiibe 12.— Load distribution along center line of rectangular wing. 

AERODYNAMIC CHARACTERISTICS 

The previous section presented solutions for the loading 
on triangular and rectangular wings flying at subsonic and 
supersonic speeds. This section will be devoted to the con- 
version of these loadings to expressions for lift and center 
of pressure. 

LEFT 

By definition the lift coefficient can lie written 



(37) 


and this will be evaluated for the various cases for which 
the loading coefficient has been obtained. 

Supersonic triangular wing. — Since the exact linearized 
value for the loading on the triangular wing flying at super- 
sonic speeds has been derived, the lift coefficient ean be 
written in the form 


a A 2 E 


(38) 


where A is the aspect ratio and E is the ell iptic integral of 
the second kind with modulus k= Vl —p i m i . 

. Subsonic triangular wing. — In the case of the subsonic 
triangular wing, equation (37) becomes 


1 C e<J d x f rai 

mc Q 2 ) Q J-mx-firfix* — y 


i 


h 



and this becomes (since A—im) 


Cl 

aA 


xofi(x 0 )dx 0 


(39) 


The numerical evaluation of equation (39) is not difficult 
since x 0 f vanishes at Xo=0. 

Supersonic rectangular wing.- -For values of Co<2*/S the 
exact value of the lift coefficient on a rectangular wing flying 
at supersonic speeds has been obtained and can be written 
in the form for jSA>l 


C L 4pA-2 
<xA~ eA ) 2 


(40) 


When PA<^1 equation (37) must be used in connection with 
equation (24) and there results 
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which reduces to 


^-£J> f'J f /l (^) is 


and this can be written in the form 


<\ p» , / 2x \ , 

ari j3ri.CoJo ** \CoPAj aX 


which becomes, if j- ; = ~-r> 
CopA 

for (3A<1 


a-u>** 


(41) 


Subsonic rectangular wing. — The equation for the loading 
on a subsonic rectangular wing, equation (14), placed in the 
formula for lift coefficient yields 



which becomes 

%=! f h(xa)dx a (42) 

<*^T X Jo 

The evaluation of equation (42) by numerical means re- 
quires special consideration since /»(*<>) approaches infinity 
at the leading edge as shown in figure 11. To this end, re- 
write equation (42 ' in the form 



(46) 


Equation (46) will be applied to the various loadings which 
have been'studied. 

Supersonic triangular wing. — The exact linearized value 
for the pitching-moment coefficient on a triangular wing 
flying at supersonic speeds has been derived elsewhere and 
can be written in the form 


On 

ocA 


3 E 


(47) 


Subsonic triangular wing. — The derivation of the pit ching 
moment on a subsonic triangular wing proceeds in the same 
manner as the derivation of lift and there results 


~~[ = — ir (48) 

This expression can be easily integrated numerically. 

Supersonic rectangular wing. — For values of pA greater 
than 1 the pitching-moment coefficient on a rectangular wing 
is given by the equation, for j8ri.>l 


C„ 6 pA — 4 

ceA 3(jS.ri) J 


(49) 


When the solution to the integral equation must be 

used and the final expression can be written 

T 

(50) 


| o /shrijtfzo +|J t M*o) dx o (43) 

and equation (B4) in the appendix (for the special case in 
which £i)=l) in the form 

1 J, f/I -> [ r+ i-a-^;j] dl °+ 

(44) 


Subsonic rectangular wing. — The equation for the pitching- 
moment coefficient on a subsonic rectangular wing follows 
in the same manner as did that for the lift coefficient. 
Hence, 

x a fi(x 0 )dx 0 (51) 

and, since the variation of x 0 fi{ x o) is as indicated in figure 13 
the numerical integration of equation (51) is simple. 

CENTER OF PRESSURE 


An application of the mean-value theorem yields 

J' f 

-r TTAETl {*- J>*» [ T Vi=b] dr °\ (45) 
L T+ RI=aU 

» 

where E§ has the modulus ki which equals p A/ -y^fl-fl 2 ) + (3 A) 1 
and where O<0<e. The combination of equations (43) 
and (45) yields an expression for the lift coefficient involving 
only the load distribution from a distance e/c a back of the 
leading edge to the trailing edge. 

PITCHING MOMENT 7 

By definition the pitching-moment coefficient about the 
apex or leading edge and based on the root chord can be 
written 


Since the pitching moment is based on the root chord, the 
center of pressure of all wing plan forms can be written 


%c.p. Efjn 

Co Ct 


(52) 


DISCUSSION OF RESULTS 

Figures 14 and 15 show the variation of the lift coefficient 
and center of pressure on triangular and rectangualr wings 
for values of P A between zero and 2. For the triangular wing, 
the differences between the subsonic and supersonic cases 
are not large in this interval of reduced aspect ratio; the 
subsonic wing develops somewhat less lift and its center of 
pressure moves forward as jS A increases. The characteristics 
of the rectangular wing, however, show a large variation in 
passing through the speed of sound. 

The subsonic rectangular wing has a variation of C L jAa 
with Ap which is almost identical with that for the subsonic 
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Fiuiihk 18.— Variation ofxo/i with xs. 



(a) Lift. 

(b) Conter of pressure. 

Fioob* 14.— Aerodynamic eharacterlstles of a triangular wing having a low value of PA. 

triangular wing. Unlike tlie triangular wing, however, the 
curve for x c .p.(c 0 on the rectangular wing shows this lift to be 
carried farther and farther forward with decreasing pA, from 
the quarter-chord position at pA= » all the way to the 
leading edge at 0/1=0. 

As the speed is further increased and the rectangular wing 
enters the supersonic speed range, the magnitude of the lift 
begins to oscillate with increasing amplitude. This continues 
until the reduced aspect ratio rises to one, after which, as 




(a) Lift. 

(b) Center of pressure. 

FimiKE 15.— Aorodsmamic characteristics of rectangular wing having a lour value of fiA. 


PA increases still farther, the value of CJAa falls uniformly 
to zero according to the expression 

xNaO-w)' (53) 

which is the exact equation for the lift coefficient given by 
linearized lifting-surface theory. The difference between 
the values of lift coefficient given by equations (41) and 
(53), represented in figure 15 by the solid and dashed lines, 
respectively, has already been discussed in the section on 
loading; the approximate solution is based on the assumption 
that the span loading at each chord station is elliptical and 
such an assumption is increasingly unrealistic for increasing 
P'A. The dotted line shown in the figure appears to be a 
reasonable interpolation between 0/1=1, the lower limit to 
which equation (53) applies, and a point where the approxi- 
mate solution given by equation (41) can be considered 
trustworthy. 

The variation of the center of pressure on a supersonic 
rectangular wing indicates that the wing is unstable for all 
positions of the pjvot point behind the leading edge for 
values of pA around 0.4, the center of pressure, in such a 
range, having moved forward of the wing leading edge. 
As pA increases past the value of 0.5, however, the center 
of pressure moves back along the wing and rapidly ap- 
proaches the midchord point, its location for a wing of 
infinite aspect ratio. 
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COMPARISON OF RESULTS FOR SUBSONIC TRIANGULAR 
WINGS WITH OTHER THEORIES AND SOME EXPERIMENTS 

DISCUSSION OP THE THEORETICAL RESULTS 

Several other published theories can be used to calculate 
the forces and moments on low-aspect-ratio triangular wings 
flying at subsonic speeds. A comparison between values of 
Ct, and center of pressure given by those theories and by 
the method of this report is summarized in figure 16. 
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Flops* 14.— Various theoretical and experimental results for triangular sings showing the 
variation of SCl. and x..,_/a with reduced aspect ratio, flA. 


The theories of Falkner (reference 13) and Weissinger 
(reference 14) are well known and will not be discussed in 
detail here. The results obtained from Weissinger’s method 
were presented by De^ Young and Harper in reference 15 
for aspect ratios equal to 1.5, 2.5, and 3.5; those obtained 
from Falkner's method were presented by Berndt in reference 
16 for aspect 'ratios equal to 1.0, 1.67, and 2.5, and by 
Falkner in reference 13 for an aspect ratio equal to 4. Both 
Berndt and Falkner used 126 vortices and 6 control points 
so the accuracy of their calculations is similar. It should 
be mentioned, however, that in order to make Berndt’s 
results consistent with Falkner’s, the values of C La given 
by Berndt have been multiplied by a factor 7 suggested and 
used by Falkner in reference 13. 

The results shown for Lawrence’s theory were presented 
in reference 17 and are based on a method originally given 
in reference 18. 

The differences between values obtained from each of the 
theories can be attributed to differences in the various simpli- 
fying assumptions used. The Weissinger and Lawrence 
theories agree on values of C Lk but disagree as to the position 
of the center of pressure. Falkner’s method and the method 
of this report yield results in good agreement for both Cl 
and eenter-of-pressure location. However, for fiA=2, the 
Weissinger-Lawrence value of C u<t is about 9 percent below 
and the center-of-pressure position about 3 percent (based 
on the root chord) farther forward than similar values 
obtained in this report 8 or by Falkner’s method. 

DISCUSSION OF THE EXPERIMENTAL DATA 

Experimental data for low-aspect-ratio triangular wings 
are given in references 13 and 19 through 22. The sections, 
aspect-ratio range, Reynolds numbers (based on the mean 
aerodynamic chord), and section thickness ratios are given 
in table ILL 

The source of the experimental values for the A =4 wing 
(given in reference 13) is some unpublished British wind- 
tunnel data. The Reynolds number is given as “high” and 
the section is not specified. Reference 21 presents the results 
of experiments made in the Langley free-flight tunnel on 
some flat-plate (% inch thick with rounded leading edge) 
models having beveled trailing edges. Hence, the values for 
the thickness ratios listed for these tests are effective thick- 
ness ratios, being, in fact, the plate thickness divided by the 
root mean chord of the wing, (2/3 )c 0 - 

The experimental data presented in table IV are portions 
of the data given in references 19 and 20 and the values 
shown in figure 16 were obtained from curves constructed 
by means of these numbers. Reference 21 presents expend 
mental results in graphical form and the slopes used in 
figure 16 were read from these graphs. It should be men- 
tioned that reference 21 also gives results for an aspect 
ratio 0.5 wing but the data presented are not sufficient to 

t a ivnrriln y to Falkner (reference 13): “It baa not been possible to establish the factor for 
all cases, bat' figures derived from a delta and other wfngs suggest that the factor can be taken 
tentatively as independent of aspect ratio, and to vary as 1+0.028 (tangent of sweepback of 
quarter chord)." 

1 It can he shown, however, that the Integral equation used by Lawrence to obtain the 
curves shown In figure IS and the one used fn this report both approach the same values of 
Ci m and r,.. Jet aa the angle of sweep goes to either r/2 or zero. In the former case they approach 
the (amUIar Jones' tow-aspectcatto+heory results, and In the latter case they yield 3r for C Lm 
and 1/3 forr«.,Vet. 
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fix the slopes of the lift and moment curves near zero angle 
of attack. Hence, no values for this aspect ratio are pre- 
sented. Finally, the values taken from reference 22 were 
read from graphs presented therein and based on integrated 
pressure distributions (measurements taken along five span 
stations). 

COMPARISON OP EXPERIMENT WITH THEORT -LIFT-CURVE SLOPE 


The experimental data shown in figure 16 (a) do not seem 
at first glance to favor either group of theories. However, 
since all these theories are based on the assumption of zero 
thickness, it is pertinent to examine the data on the basis 
of thickness ratio. A correlation on this basis is presented 
in figure 17. It shows that with decreasing thickness ratio, 
the ratio of the experimental values to those predicted by 
the method of this report (in the range 0 <0A<2) or by 
Falkner’s method (in the range 2 <(SA<4) tends to unity. 
Of course, tills trend is not conclusively bomo out by these 
comparisons and should be tested by more experimental 
measurements in the lower thickness ratio range. 

It is interesting to notice that a good approximation to 
the results calculated in this report and by Falkner’s method 
is given by the equation 


ft. 


2irri 
pA + 2 


(54) 



Pini’RZ 17,— Variation of (CAJ „,/( Ci.) u.' the ratio of the measured lift -curve slope to that 
given by the theoretical methods of Falkncr or this report with thickness ratio. 


where 


P= 


wing seinipeiimeter 
wing span 


Equation (54) was derived by It. T. Jones (see refercnco 23) 
as a first-order correction to the value of C La given by 
lifting-line theory for wings having elliptic plan forms. It 
has been found, however, to proride a good estimate for the 
aspect-ratio correction to wings of various plan forms.® For 
the particular case of the triangular wing in compressible 
flow, equation (54) becomes 


PC L - 




■+¥+ 


V 


>+(¥)’ 


(55) 


and its variation for 0<£A<4 is shown in figure 16 (a). 

COMPARISON OF EXPERIMENT WITH TUEORV-CENTER-OF-PRESSURE 

LOCATION 

The comparison between experimental and theoretical 
values of center-of-pressure location is shown in figure 16 (b). 
In weighting the experimental points shown, it should be 
remembered that the tests recorded in reference 21 were 
made on a flat plate with a beveled trailing edge. When 
such a wing is at a positive angle of attack, the bevel causes 
the flow to separate at the surface discontinuity producing 
an effective upwardly deflected flap. The result is a loss of 
lift on the rearward portion of the wing and a forward shift 
in center of pressure. Hence, the experimental values for 
Xc.pje 0 taken from reference 21 should bo low relative to the 
values given in the other references. With this taken into 
consideration, the experimental values shown in figure 16 (b) 
are in fairly good agreement and again favor, at least, the 
theoretical results based on Falkner's method. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Nod. 28, 1952. 


* This approximation is suggested In an article by R. T. Jones and Doris Cohen, “Aerody- 
namics of Wings at High Speed,” to be published by the Princeton University Tress In 
section A of the book entitled Applied High-Speed Aerodynamics, volume VI of the High- 
Speed Aerodynamics and Jet Propulsion Series. 



APPENDIX A 

EVALUATION OF SPECIAL INTEGRALS 


THE INTEGRAL Z, 

The evaluation of I t will be discussed first for the case in 
which x^>Xi and second for the case in which x<C*i- 

CASE I, x>xi 

It is possible to write I t in the form 

f'jg2S=i, cad 

where Mo=nUx— and fh=m{x+xi). The linear term in 
the lower radical of the integrand can be eliminated by the 
transformation »j = (<r+5t)/(l+t) > and the integral becomes 


j j-»_ _ . r- 

VG*o“ <r) (<t — JZZi 


vo 

PQ— ^ 


where 


V(M, 1 +<rVm»)+G x 0 »+3’gW dt 

a+ofr+^i 

r ~ Mo(mq— /3*m VQ +mo V(l + P*™*) (Mo^+ft^VJ) 


and 


-i^oOio— j3 2 mVi)— 


(A2) 

(A3) 

,(A4) 


The expression for <r and S may be combined to give the 
useful identities 

Mo I= — <r5j3 s m* 

( mo — a-) ( J — Mi ) + (mi — a ) (5 — Mo) = 0 

Using fundamental properties of even and odd functions, 
equation (A2) may be reduced to the form 

T-o I rf+a'pm* fT/- - 1 tfr \ 1 

_ \ Bk r <rk f J 


/t ,2 +tv 

V-T=3- 


bv the substitution 


du 


<■>= 


Mi~ 5 


Mi— 0" 


and where 


M_ (mi— 0-)^ 
(ml 2 — 0-5) (5— <r) 


P 5 / Mi' — g V 

£'* ff\Mi — 5/ 


By introducing the Jacobian elliptic functions in the 
transformations w=cnu, the integral reduces to 


r = 2 /_M o 

1 (m — <r)(<r — 


where 


f*( n a a y ** ) 

— cr)(<r — mo) Jo \l+i'i*aa*u 1 +v 1 a in t u/ 




k fl +~ F 
0 


>0 


dnhidu 
(A5) 

CA6) 


and 


r 2 


- ~u 


->o 


XA7) 


F*+ " F 

O' 

The integration n ay now be completed and 

T _2 / mw^™* r i ^TF Vi \ 

1 &A (mi — o-)(ff— Mo) L A \yt>*-fP/ 




Mg I 

f —Tf— ^”1 

i Mi 1 

W**+*vJ 


where 

G(;x)=EF'(x ) +i£E'(x)-iiy / (x) 
the modulus of the elliptic integrals being £■ or F. 

CASE 2, x<x. 


(AS) 

(AS) 


The procedure for obtaining the solution for h in tliis case 
is identical to that followed in case 1 except that in order to 
fulfill the contition that £> — 1, <r and 5 must be defined in 
the following manner 


_ — mq(mq— ^ m^O+MoVtl+^^CMo^+ff^Vi*) / yo) 

j3*m’(Mid Mo) 


— Ma(Mo— ^mVi)— M oVII+iS^^Cmo^IS^Vi 2 ) 

FmKjii ~F Mo) ^ 

In this case it can be shown that in equation (A5) V])>0 
and Vi<— 1, and the solution for h is 
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/,=2 

em VSKwWw)> 

i u^3±^i k \ : 

8 V {Ml — <r)(tr— Mo) ; 
where 0(x) is defined as in equation (A9). 

THE INTEGRAL 7 2 
Writing I 2 in the form 

r_o f' /(s-Xi) s +yi*i3*_,.. 

V~ ^=7.*“ _%l 


(A 12) 


(A13) 


and setting 


—=cnu, k = — == ^ = 
s V(®— a?i)*+^V 


7 S can be integrated to give 


* g dntudU—ZpsE 


I 2 =2ps j 

THE INTEGRAL / s 

CA 

Writing J s in the form 


ok k 


CASE I, 


(A 14) 


V Mo 2 2 
0 q 2 ” 

(mi-vXh-wT 


and making the transformation 

2 , y — mo 

Tl -W I Mo 
iTr 
“0 


■s'?™ 




where 


fc». 


20o Ml— Mo 


1—0 O Mi0o + Mo 

reduces equation (A15) to the form 


(A 15) 
(A 16) 
(A 17) 


! k*sn 2 u 


/ 3 -(l +0o)^ ^- mo) /"X x [ l-f 

\ 20 o 
6 0 du \ 
l_i^op S n s a J 

The integration may now be completed and 

J (1 _ e 0 ]K^ 1 - - ^ 1 +^- ■ 

\ Mi— Mo J 


00-y/l"- 


1 0Q jr,j 


20 o 


(A18 > 


CASE A 0<e»<5t^i 


In this case 


/j — 0o I -tl d* dij 

*° V (Mi 7) (v — Mo) V 


(A 19) 


and the transformation 

S7l*U= 


is made where 


2. V-Mo 
1 — 0o 1 Mo 

,+ 5; 


Mo. 


l-0o A<1 + 0 


2 Mi— Mo 

Equation (A19) then becomes 

iz={i+6 0 )J I( - 2 ^- J f Jf — pr f*— iH r — ^ 

V 0 o(mi-Mo)^Jo l+L^sn’u J° 

The integration can be completed so that 

/,=<1+W V5^{^(l-«.+2«ii i-1 ) + 

/ 20p(l — 0Q 8 ) ' ° (V 1-00+200^) . ,, 

V 1 — 0o+20ok 2 l4~0o 

0 O / 2(l-0„ 2 ) /l /l^V 

1 -f 0 O V 2£ l - 1 + 00 L A E u V 2 / 

(W¥)| 




THE INTEGRAL /, 

CASE 1, 0<Xi^JT-Jfc 

The integral h can be written 


(A20) 


< A2 » 

When the transformation snu—yi/s is made, the expression 
may be integrated to give 

f i =2(x~xi)j' dn ! udu=2(x—Xi)Ei (A221 


where 


*!* 


8 2 jS s 


(l-Kl) 1 
CASE 2, X-Pt<X t £X 

In this case /* can be written 

I~Ti 


(A23) 


/ - 2 X < <**> 

The transformation snu=^y l j{x—x\) applied to equation 
(A24) yields 




where 


and 


— ^ ■ cn t udu=2k i (x— X\)B% 

(A25) 

*2 

pa 

(A26) 

D Ei—kt^Kt 

Bi ~ — s? — 

-(A27) 
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Since the integral h is a function only of the ratio xifx, 
equation (10) can be written for w=v.' 0 


‘-if . 1 (ED 

where £=xi/x. It is novr assumed that fi (x£lc 0 ) may be con- 
sidered constant over small intervals. This reduces the 
solution of the integral equation to the elementary problem 
of solving a system of simultaneous algebraic equations. 
On the basis of such an assumption, equation (Bl) becomes 

a 


■a -1 


where 


So(©= 




i-e 


m; 0<|< 


Soi5) = Yy^ I'd©; 


<£<^ 


(B3) 


The function g 0 (© was calculated by numerical integra- 
tion and systems of simultaneous equations were obtained 
for values of n equal to 3,. 6, and 9. Solutions were found 
using the Gauss-Seidel method (for which the simultaneous 
equations were well suited). 

SUBSONIC RECTANGULAR WING 

Substituting the value of Jj given by' equation (A14) into 
equation (16), one has for w=uc 0 


where 

and where 

k= 




(B4) 


0s 


0A 


VCx-rd’-H/SV 




+jS i A 1 


A satisfactory numerical solution of equation (B4) requires 
the solution of the system of simultaneous equations of the 
form 


1 = 


( V -°)+ 2 [a (£> * firv) 


+ 


a'Xtv-O} <B5) 

The convergence of the solutions to equation (B5) is indicated 
in figure 18 where’the value of n was successively taken to be 
3, 6, and 9. 


1.2 

i.O 

.8 

» 

.6 

A 

.2 

O, 



\ 

\ 

\ 










A 




■ 







k 

\ 

\ 


O 3 Intervals 
* 6 Intervals 
□ 9 Intervals 






\ 

\ 

V 







/ 

\ 

Yr 







-1 



h 












IT 5 





.2 


A 


.6 


.8 


1.0 


Fir.r** 18.— Solutions obtained for the variation of fi svtth ri using S, S, and 8 Intervals in 

equation (B5). 

SUPERSONIC RECTANGULAR WING 
Equation (26) can be written when ic=ic 0 as 


where from equations (A22) and (A25) ' 


(B6) 


(2 

l!iU T 




E x 


X~Xi 

03 


>lj 


and where 




(B7J 


*i= 


S0 


X — Xi 


kt* 


X—X! 

S0 


By the application of the trapezoidal rule for numerical 
integration to equation (B6), it is possible to write fi(xfs0) 
explicitly' as 

<W> 

where A(x/j3e) is the interval of the trapezoid. 
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APPENDIX C 


EVALUATION OF S 

C '(insider tlie integral equation 

■' (C1 ) 


2rI „_ p ,, (•' 

do \Xo—Xi J 


X 0 —Xi 


In order to study /i(zj) in the neighborhood of the origin, 
first set J' 2 =x 0 f so that equation (Cl) becomes 


_ C 1 , P /Iu Ui(S)/i(*of)^ 


-X 


Vi-? s 


i-i 


(C2) 


then assume that, for small values of x?, (i-i) can be 

expressed as 




From this one has 


2 „-«» ^ r, r 43+ r ^ 

j*-h0 *o*L Jo Vl-f* Jo 1 — ^ J 


which immediately implies 

r e-dt , r (£)<** „ 

"J. 7f=? + j. “T=f 0 


(C4) 


Tlie value of 3 in the LnU'rval 0</9A<2 that satisfies 
equation (C4) was determined by numerical integration. 
The results are shown in figure 8. 
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TABLE I.— VALUES OF CONSTANTS FOR EQUATION (35) 


Const-N. 

0 

0J. 

OH 

a* 

I. 0000 

L0040 

1.0513 


0 

.3151 

.0166 

Oj 

0 

—.6422 

-.mi 

y 

0 

.010 

.022 


TABLE II.— VALUES CALCULATED FOR RIGHT SIDE OF 
EQUATION (33) USING EQUATION (35) AND TABLE I FOR 
/. (*«) 

[Exact solution would, In each cast, yield 1.0001 


n 

(drap-0.1 

(fim)'- 0A 

0.1 

0-985 

0.991 

.2 

L 000 

1.000 

.3 

L011 

LOOS 


L 022 

LOW 

.5 

L022 

LOOS 

.6 

L016 

L007 

,7 

L010 

.999 

.8 

L000 

L000 

.0 

.993 

.997 

.95 

LQ0G 

LOOO 


TABLE III.— RANGES OF PERTINENT PARAMETERS IN 
EXPERIMENTS 


Reference 

number 

Section shape 

Thickness 

ratio 

(percent) 

Reynolds 

number 

Aspect ratios 

13 

to 

10 

« 

4 

19 

Swedish 

10 

3X10* . 

MM3 

20 

FFA 104-5106 
N'ACA 0012 

12 

L7 to 3X10* 

L4/3A3 

21 

(') 

(2 to 4)1 

.4 to 9X10* 

22 

MAC A 0012 

12 

2.4X101 


TABLE IV.— EXPERIMENTAL DATA ON LOW-ASPECT-RATIO 
TRIANGULAR WINGS TAKEN FROM REPORTS BY 
BERNDT (REFERENCE 18) AND LANGE AND WACKE 
(REFERENCE 20) (SEE TABLE III) 


(a) Reference IS (C. referred to id end taken about point {aback from apex) 


A- S/2 

A-S& 

i 

1 

a 

Cl 

D 

a 

Cl 

B 

at 

Ct 

Cr. 

-3.7 

-0.108 ' 

0.029 

-AS 

-a 165 


-A 2 

-a 097 

0.022 

-Z7 

-.122 

.022 

-3.8 

-.130 


-3.2 

— .072 

.018 

— L 7 

—.073 

.014 

-2.8 

-.096 

OBtiTvB 

— 2.2 

-.049 

.012 

-.6 

-.027 

.006 

-L7 

-.060 

am 

— 1.2 

-.028 

.007 

.4 

.020 

-.002 

-.7 


.008 

—.2 

-.005 

.002 

L5 


-.009 

.3 



.8 

.018 

-.003 

2.5 

.115 

-.017 

L3 

.043 

-.005 

1*8 

.042 

-.009 

3.5 


-.024 

2.3 

.078 

-.012 

2.8 

.064 

-.014 

4.6 

.205 

—.031 

3.4 

.114 


3.8 

.068 

-.030 




4.4 

.149 


4.8 

.110 

—.025 


(b) Reference 20 (C. referred to {a and taken about point |C| back from apex) 


A-l 

A-VZ 

a 

Cl 

C» 

a 

Cl 

c- 

—5.84 

0 

5S4 

-0.1258 

-.0002 

.1255 

-0.0047 
--OO0O (sic) 
.0067 

-5.81 

0 

5.80 

— (L 152 
.002 
.156 

-0.0101 

0 

.0089 

A-2 

A-Z 

CL 

Cl 

c- 

CL 

Cl 

Ml 

-5.71 

0 

6.76 

-0.2375 

-.0180 

.1961 

-a 0177 
-.0007 
.0180 

-5165 

.02 

5.09 

-a 2830 
-.0162 
.2561 

-a 0285 
.0002 
.0311 


' See text tor explanation or qualification. 





































